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A qualitative analysis of a scalar-tensor cosmological model, with an exponential 
potential for the scalar field, is performed. The phase diagram for the flat case 
is constructed. It is shown that solutions with an initial and final inflationary 
behaviour appear. The conditions for which the scenario favored by supernova type 
■ la observations becomes an attractor in the space of the solutions are established. 
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O ! 1 Introduction 

o 

The results from the high redshift supernova type la observations |2| indicate that the 
7 5 " 1 , Universe is in an accelerated expansion regime. This means that the matter content of 

the Universe must be dominated by an exotic fluid whose pressure is negative with p < 
— |. Generally, the cosmological constant, with p = — p, is assumed as the most natural 
candidate to represent this exotic fluid. However, the theoretical problems concerning a 
cosmological constant, somehow connected with its interpretation as a vacuum energy ||, 
make other possibilities very attractive also. 

Some kinds of topological defects can also lead to an effective equation of state in 
the searched range. Non-relativistic domain walls, for example, implies p = — |p. But, 
a more fashion theoretical proposal is the so-called quintessence model [Q, which is a 
minimal coupled scalar field with a slow evolving potential. It has been argued that such 
kind of potentials may originate from supergravity models, and the accelerated solutions 
correponds to an attractor || . These kind of models may be tracked back to the reference 
H, where they have been proposed in order to solve the dark matter puzzle. Quintessence 
model makes use of a large variety of forms for the potential: exponential^], polynomial 
combined with exponenential J5], hyperbolic sine or cosine @ or even a double exponential 
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In a recent work ||, it has been showed that an exponential potential model admits as 
particular solution the typical perfect fluid solution of the Friedmann-Robertson- Walker 
model, for any value of the barotropic equation of state parameter a (p = ap). Hence, 
inflationary power-law solutions are covered by this model. One of our interest here 
is to study the status of this particular solutions, verifying to which extent they are 
attractors, and if it is possible to obtain models where the expansion is initially non- 
accelerated, becoming accelerated later. Due to the fact that the standard cosmological 
model needs an inflationary phase in its primordial phase, it would be also interesting 
to have models where an initial and final inflationary phase occur, with an intermediate 
non-inflationary behaviour, in order not to spoil the nucleossynthesis achievements and 
structure formation. 

In order to do so, we will perform a qualitative analysis of this model, first in the scalar- 
tensor model, and secondly with the scalar-tensor model coupled to ordinary matter. Some 
studies have already been made for these cases in the literature [|K], O, |I2[ . However, 
in |10| it has been mainly verified that the flat case is an attractor for a great variety of 



exponential factor, while in ||TT], O the analysis, with respect to this exponential model, 



was mainly dedicated to the identification of the nature of critical points. 

In the present case, we will be interested in mapping completly the solutions for k = 0, 
since it has already been shown that k = is an attractor with respect to k = ±1. Another 
reason is that k = seems to be favoured by observations [^] . The phase diagram will be 
constructed and the positive energy and inflationary regions will be identified. It comes 
out that the both kind of desirable models, as described before, appear and the power-law 
particular solutions are indeed attractors of the physical acceptable solutions. 

In next section we describe with some detail the scalar-tensor model with exponential 
potential. In section 3, the phase diagram is constructed. The coupling to ordinary matter 
is discussed in section 4, while in section 5 the conclusions are presented. 



2 The scalar-tensor model 

A minimal coupling between gravity and a self-interacting scalar field is represented by 
the lagrangian, 

L = V=ff[i2-M ; " + 2F(0) 
The field equations are 

Rnu - ^g^R = <P^<P;u - ^g^^ + g^V(4>) , (2) 
□0 = -V\<j>) , (3) 

where the prime means derivative with respect to (p. Inserting the FRW metric 

r dv^ 

ds 2 = dt 2 -a 2 {t) — + r 2 {d8 2 + sin 2 6d<j) 2 ) (4) 

. J. KV 
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where k — 0, —1, 1 corresponds to a flat, open and closed Universe respectively, it results 
the following equations of motion: 



3(-) 2 +3A = If+VM , (5) 

-2^-(^) 2 -A = i^-y( ) , ( 6 ) 

a \a/ a 2 2 

+ 3-0 = -V'fo) . (7) 
a 

These equations are connected by the Bianchi identities. If the potential term is written 

as 

V(<P) = l-^\ 2 exp(TVW+aj<P) , (8) 
6 (1 + ay 

where — 1 < a < 1, the equations of motion (|^,|) admit, for the flat case, the solutions 

a{t) = a tm , 0(t) = ±— ===]nt . (9) 

V3 1 + a) 



Hence, the potential (H) leads to the usual perfect fluid solutions of Einstein equation 
as a particular case. When a > — ~, these particular solutions represent an expanding, 
decelarating Universe; if a < — ~, the scale factor describes an expanding, accelerating 
(inflationary) Universe. The limiting case a — — | corresponds to a oc t, that is, a > 
and a = 0. So, the perfect fluid solutions may be mimitized by a scalar-tensor model with 
a suitable potential. In it was shown that this scalar-tensor model allows to get rid 
of instabilities that appears in the perfect fluid models, at perturbative level and in the 
small wavelength limit, when a < — |. 

However, it is important to know if these particular solutions represent an attractor of 
the space of all possible solutions of this model. Moreover, the solutions @ describe an 
Universe that is always inflationary (a < —1/3) or always non-inflationary (a > —1/3). In 
a realistic Universe, the behaviour of the scale factor must change with time. If we accept 
the inflationary paradigm for the early Universe and if the results for the value of the 
decelerating parameter q = — ^ obtained through the measurement of the supernova type 
la are confirmed, then the most realistic model should have an initial and final inflationary 
behaviour; between these two inflationary stages, a non-inflarionary behaviour must take 
place, allowing the formation of light elements, through the primordial nucleosynthesis, 
and formation of local structure through the gravitational instability mechanism. 

The model described by ([[]), with the potential given by (||), is rich enough to allow 
solutions more general than (|9|). In order to exploit all richness of this model, verifying if 
realistic models as described before are possible, we will perform a qualitative analysis of 
(^.[|). Compactifying the phase space of all possible solutions on the Poincare's sphere, 
and delimiting the regions corresponding to the positivity of the potential V((f>), as well 
as the inflationary type solutions, we will indentify the main features of the cosmological 
models resulting from ([||||7p. 
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3 The phase space diagram 



In order to perform a qualitative analysis of the system described above, we define the 
new variables 

a 

x = 4> , y = - , (10) 

a 

and we set B = ±a/3(1 + a) and V(<p) = V e B( ^. Hence, equations (f||6|,|7]) lead to the 
planar, homogenous, autonomous two dimensional system 

x = —x 2 — 3xy — 3By 2 , (11) 



2 
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k , (12) 



subjected to the condition 



V{(j>) = y 2 - \x 2 > . (13) 
o 



This system admits in the finite region of the plane(x, y) an unique degenerate critical 
point x — 0, y — 0. It corresponds to the Minkowski space. Its eigenvalues are all zero. 

Now, the invariant rays in this plane are characterized by the solutions x = Xy. 
Inserting this relation in ( |TT| , |r2"D it results a third order polynomial relation for A: 

X 3 + BX 2 -6X-QB = <=^> (X + B)(X 2 -6) = <=> X 1 = -B , A± = ±VQ .(14) 

Hence, there are three invariant rays for B ^ ±v6. Two of them are independent of the 
value of B while the third one depends on B and, consequently, on a. This third inviariant 
ray correspond to the solution @. Let [XY], [CC] and [AA'] denote respectively the rays 
x = —By, x = —VGy and x = V~Gy. 

Inserting the expressions for these invariant rays into flT2p the direction of evolution 
along them can be determined. In order to complete the analysis, the critical point at 
infinity must be found and their nature studied [fT4"| . To do this, we project the plan (x, y) 



into the Poincare's sphere, introducing new variables z,u,v, such that x — — , y — — 

z z 

subjected to the condition u 2 + v 2 + z 2 = 1. The new system reads: 

— = bz — cv , (15) 
dr 

— = cu — az , (16) 
dr 

dz , 

— = av — ou , (17) 
dr 

where a = \u 2 z, b = z(^u 2 — 3uv — 3Bv 2 ) and c = — 1-u 3 — ^u 2 v + 3uv 2 + 3Bv 3 . The 



new parameter r is defined by the equation (JlTf) - The points at infinity are obtained for 
z = which corresponds to the equator of the Poincare's sphere. The critical points are 
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determined and their nature (repulsive, attractive or saddle points) is characterized in the 
usual way [[14]. The coordinates of these points are 



A(u,v,z) = \\~,— ,0 , A'(u,v,z) 





X(u,v,z)—\—= ,—z , , Y(u,v, z) 

The final phase diagram contain six regions separated by the three invariant rays. 
Three of the six regions may be obtained from the other three by inverting the time: 
t — > —t. In order to interprete the behaviour and the nature of the curves, we complement 
this diagram with some other physical considerations. 

The complete diagram depends on the value of B. However, before establish these dia- 
grams two physical requirements will be introduced. The first one concerns the positivity 
of the contribution of the energy of the potential term. The potential V{4>) is negative in 
the region interior to the rays such that A = In this region the variable x can never 

be zero, since in this case the variable y becomes imaginary. 

We will be also interested in identifying the regions where inflation can occurs. This 
implies to require a > 0. Hence, y + y 2 > 0. Due to (|12D , this implies — \x 2 + y 2 > 0. 



Hence, the inflationary region is bounded by the rays y = ±^x. These rays will be 
represented by [PP'] and [QQ']. The inflationary regime occurs in the region bounded by 
these lines and containing the y axis. 

Hence, in these phase diagram the two invariant rays (from three), the condition of 
positivity of the energy, and the condition to have inflation are fixed and independent of 
B. Only one invariant ray depends on B; it corresponds to the particular solution (§). 
When a = 1, this invariant ray coincides with one of other two. 

The complete diagram is displayed in figure 1 for — y/2 < B < 0. As B changes 
from —a/6 to VQ, the invariant ray [XY] moves from [CC] to [AA']. When a = — 1, the 
invariant ray [XY] coincides with the y axis. 

If — | < a < 1, (— y/6 < B < —\pl or y/2 < B < \/6), the particular perfect fluid- type 
solutions are non-inflationary. There are two physically interesting kind of trajectories 
in this case. Those connecting the critical point at infinity C to the origin, approaching 
asymptotically the invariant ray [YX], start with a free scalar field behaviour (a oc t 1 ^ 3 ), 
and coincide later with the corresponding perfect fluid solution; these curves describe 
non-inflationary Universe. However, the curves connecting the critical point at infinity 
Y to the origin, start with a non-inflationary perfect fluid behaviour, becomes latter 
inflationary and assymptotically coincides again with the corresponding non-inflationary 
perfect fluid solution. 

When — 1 < a < — |, which corresponds to — y/2 < B < y/2, the perfect fluid- type 
particular solutions corresponding to the invariant ray [XY] are inflationary. Again, we 
have two kind of interesting trajectories. The solutions starting from the critical point 
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at infinity C begins again with a free scalar field, non- inflationary, behaviour and then 
evolves towards an inflationary regime as they approach the origin. Those solutions 
starting from the critical point Y, have initially an inflationary behaviour, becomes later 
non-inflationary and again become inflationary as they approach the origin along the 
invariant ray [XF], There are also solutions that begin at Y, ending in the origin, 
remaining always inflationary. 

These last trajectories are the most interesting one since the standard cosmological 
model requires an initial inflationary regime, as well as a final inflationary regime, if the 
supernova type la measurements are confirmed. 

When B — (a = — 1) all y axis becomes singular and corresponds to the de Sitter 

particular solution with different values for the Hubble factor — . All physical acceptable 

a 

solutions go to or come from one point of this singular axis. It is important to notice 
that, for any value of B, the invariant ray [XY] is an attractor. 



4 Coupling the scalar field to ordinary matter 

The inclusion of ordinary perfect fluid matter, with a barotropic equation, can be made 
in a quite direct way. In this case, the equations of motion read 

~) 2 = 8nGp + \tf + V{<j>) , (18) 

- 2 --(-) 2 = SnG^p + U'-V^) , (19) 
a \a/ 2 

+ 3-0 = -V^<P) , (20) 
a 

p + 3(l + 7 )-p = . (21) 
a 

In these expressions, the barotropic equation of state p = 7p was explicitly employed. 

The equations of motion above may be recast in the form of a three-dimensional 
dynamic system: 

x = Bz + —x 2 - 3xy - 3By 2 , (22) 
y = - l -^z + x 2 , (23) 

i = -3(\ +1 )y Z . (24) 

The definition of x and y are the same as before and z = 8irGp. For z = 0, we go back to 
the two-dimensional system (pilIIS). 

The study of a three-dimensional dynamical system is more involved. However, we 
may obtain its general features following closely the procedure employed in the two- 
dimensional case. First of all, we remark that the system (p2|p3|,p^) has one critical 
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Figure 1: This diagram gives the Evolution of different solutions throught inflationary or non inflation- 
ary regime. It shows the reheating phenomena for some solutions in the sector OCX and YOA 
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point, as before, represented by x = 0, y = and z = 0. It corresponds to the Minkowski 
space. But, in order to obtain the invariant rays and the critical points at infinity, we 
must perform some suitable "cuts" in the three-dimensional phase diagram, reducing it 
to an ensemble of two-dimensional system which can then be completely analyzed. These 
"cuts" correspond to projections of the three-dimensional trajectories on some planes. 

The most simple and natural projections correspond to impose x = 0, y = or z = 0, 
respectively. The last one leads to the two-dimensional system analyzed in the previous 
section. On the other hand, when we impose x = (no scalar field) or y = (no gravity), 
the resulting system can be completly solved. Specifically, these hypothesis lead to 

x = -> y = J| -> a oc t^rW ; (25) 
y = -> x = z = (a > -1) . (26) 

The first case corresponds to the gravity perfect fluid system, whose solution is well known. 
The second one corresponds to the trivial case: the Minkowski solution can exist only if 
matter and the scalar field are absent also. 

The fact that those two cases are completly solved (one of them through the trivial 
solution) just mean that if the solution is initially in a plane x = or y = 0, it remains 
there. Hence, the most interesting solutions, with non trivial solutions for x, y and z, are 
those outside these planes: however, their projection on the z = behaves as described 
in the previous section, and all analysis performed before remains. 



5 Conclusions 

A self-interacting scalar field coupled to gravity is considered as a good candidate to 
describe the dominant matter content of the Universe today, leading to an accelerate 
expansion. The potential term in general is taken as an exponential function of the scalar 
field, or a combination of polynomial and exponential functions. Here we have exploited 
a model where the potential is just an exponential function of the scalar field. In || it 
was shown that such potential term may lead to power-law solutions for the scale factor 
typical of a perfect fluid gravity system, with an arbitrary barotropic equation of state. 
Hence, inflationary power-law solutions are included in this scalar-tensor model. 

In the present work, we performed a dynamical analysis of that model for a flat 
Universe. It was verified that the power-law particular solutions act as attractor in the 
space of the allowble solutions. Moreover, it was identified, for some exponential factors, 
solutions with an initial and final inflationary behaviour as well as solutions with just a 
final inflationary behaviour. They coincide asymptotically with the particular power-law 
solutions. 

The restriction to k = was made due to two reasons mainly: in , it was shown that 



in an exponential potential scalar-tensor model, the flat case is in general an attractor. 
Moreover, the CMB anysotropy observations favor a flat Universe [jnj. The goal of the 



present work was to present a complete phase diagram description for this particular, but 
very important case. 
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The results showed the richness of the model, and that such a simple self-interacting 
scalar field may lead to scenarios that are consistent theoretically and can be in good 
agreement with observations. In particular, if we admit a potential of the type V(<f>) = 
Voe ±B ^, with B = 1, the scale factor behaves asymptotically as a oc t 2 , which is one of 
the most likely behviours for the scale factor as can be infered from the supernova type 



la observational programs p5 
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